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SUMMARY 

This  report  outlines  a technique  for  deriving  a set  of  initial 
error  states  from  an  error  covariance  matrix.  The  need  for  an 
ensemble  of  state  errors  arises  in  a Monte-Carlo  error  model  for  the 
re-entry  portion  of  flight  of  a ballistic  missile.  Typically  the 
boost  and  vacuum  coast  part  of  a trajectory  is  accurately  modeled 
with  a covariance  propagation  scheme.  During  reentry,  however, 
errors  become  nonlinear  in  time  and  are  highly  cross  correlated 
with  one  another  and  interact  strongly  with  the  dynamics  of  the 
vehicle  flight.  To  accurately  model  the  errors  during  reentry  a 
Monte-Carlo  approach  is  often  required.  The  interface  between  a 
covariance  propogation  scheme  and  a Monte-Carlo  model  requires  that 
an  ensemble  of  initial  condition  errors  be  extracted  from  an  arrival 
covariance  matrix.  The  errors  must  be  properly  correlated  and  have 
the  proper  distribution  of  magnitudes.  The  method  described  herein 
is  one  technique  for  defining  the  interface  between  the  two  types 
of  error  models. 

Since  the  original  draft  of  this  report,  the  author  has  found  that 
a method  with  equivalent  results  has  been  used  to  determine  wind 
profiles  from  a matrix  of  wind  correlation  coefficients.  (e.g. , 
Hankerson,  S.  H. , "Wind  Profiles,"  NWL-TN-G-4/72 , Feb  1972).  It  is 
likely  that  other  uses  for  the  general  method  can  be  found,  especially 
in  the  field  of  system  error  modeling. 

The  present  work  was  sponsored  by\^^  Navy's  Strategic  Systems 
Project  Office,  Mr.  Roger  Stanton,  SP-5^232,  monitor,  as  part  of  the 
Improved  Accuracy  Program  (lAP) . Progipim  guidance  of  Dr.  J.  Goeller, 
Acting  Re-Entry  Technology  Co-ordinator,'  Naval  Surface  Weapons  Center, 
is  gratefully  acknowledged. 
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I.  INTRODUCTION 


Often  it  is  desired  to  determine  the  performance  of  a 
physical  system  in  an  average  or  statistical  sense.  That  is, 
for  many  repetitions  of  a phenomena,  what  are  the  operating 
limits  within  which  a system  may  be  expected  to  perform  and  how 
are  deviations  from  nominal  behavior  distributed  about  the  average 
of  the  many  repetitions.  A study  of  the  statistical  behavior  of 
such  physical  phenomena  is  here  referred  to  as  system  error 
modeling. 

There  are  two  common  methods  of  constructing  system  error 
models.  The  first  method  is  the  deterministic  or  Monte-Carlo 
f method.  This  method  attempts  to  mathematically  model  the 

physical  laws  which  affect  a system  and  allow  the  parameters 
within  the  physical  laws  to  vary  in  a particular  manner  within 
their  limits  of  uncertainty.  The  parametric  variations  allowed 
have  a certain  randomness  associated  with  them  since  the  exact 
values  of  the  parameters  are  never  precisely  known.  Correlations 
between  the  parameter  errors  are  permitted  if  they  can  be 
determined  a priori.  The  system  is  allowed  to  operate  within 
the  framework  of  the  Math  model  for  many  realizations  of  the 
event  being  studied.  After  exercising  the  model  in  this  Monte- 
Carlo  mode,  the  system  states  for  the  many  realizations  of  an 
event  can  be  analyzed  and  statistical  variations  of  state 
amplitude  and  distribution  can  be  determined.  This  statistical 
description  of  the  system  is  the  result  being  sought. 

The  second  method  of  system  error  modeling  is  the 
covariance  propagation  scheme  or  analytic  model.  This  model 
I.  represents  the  uncertainties  or  errors  in  the  state  of  a 

' system  in  a covariance  matrix.  The  elements  of  this  matrix 

describe  the  magnitude  of  state  errors  and  their  correlations 
? with  one  another  in  a statistical  sense.  The  analytic  error 

j * model  propagates  this  covariance  matrix  with  time  to  describe 

. how  the  state  errors  change.  Since  all  of  the  desired  statistical 

; information  is  available  in  the  covariance  matrix  it  need  only 

; . be  propogated  once  for  each  event.  This  method  has  the  advantage 

of  reduced  complexity;  however,  it  is  not  as  versatile  as 
the  Monte-Carlo  method  and  determination  of  the  proper  way  of 
propagating  a covariance  matrix  is  not  always  a simple  matter. 
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In  some  instance  it  may  be  desirable  or  necessary  to  mix 
the  two  methods  of  error  modeling  for  a given  physical  system. 
If,  say,  the  initial  part  of  an  event  is  modeled  with  the 
analytic  error  model  and  the  latter  stages  with  a Monte-Carlo 
model,  then  how  are  the  two  models  interfaced  at  the  point  of 
transition?  What  technique  will  be  used  to  calculate  the 
ensemble  of  initial  deterministic  state  vectors  from  a terminal 
covariance  matrix  in  order  to  begin  the  Monte-Carlo  process? 

The  following  analysis  describes  one  method  of  defining  the 
interface. 


II.  LINEAR  REGRESSION 

A.  Two  Variable  Regression 


Given  an 

(n  X n) 

symmetric  error  coveuriance  matrix: 

^11 

^12 

®13  • 

• * ®ln 

^21 

• 

®22 

®23  • 

• * ®2n 

• 

Snl 

Sn2 

Sn3  • 

• * ^nn 

The  diagonal  terms  of  this  matrix  represent  the  variance  of 
error  in  a state  property.  If  a state  property  is  quantified 
by  a known  part  plus  an  uncertain  zero  mean  error  in  the 

state  property  say  Y^,  then  the  property  Y may  be  expressed  as; 


Y 


+ ^u 


(la) 


The  term  of  the  error  covariance  matrix  then  represents 

the  statistical  variance  of  Yu  for  many  realizations  of  the  state. 


S 


11 


(1) 


N is  the  large  number  of  realizations  of  Yu.  Similar  definitions 
exist  for  the  other  (N  - 1)  states  implicit  in  the  covariance 
matrix.  The  off-diagonal  terms  of  the  matrix  cire  the  covariances 
of  the  (N)  states  defined  by: 


(2) 
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The  magnitude  of  the  covariances  indicate  the  degree  to  which 
state  errors  are  dependent  on  one  another.  To  put  this  dependence 
in  a nondimensional  form  a correlation  coefficient  can  be  written: 


’12 


12 


^®11  ®22^ 


(3) 


In  this  nondimensional  form,  a correlation  of  +1  indicates  perfect 
correlation  between  errors  ^u^  and  ^U2,  and  a correlation  of  0 
indicates  no  correlation  exists.  This  definition  of  correlation 
coefficient  is  applicable  only  for  linear  correlation  analysis. 

A non  linear  correlation  may  exist  for  ri2  ® present 

method  assumes  only  linear  correlations  exist  among  the  state 
errors . 


If  the  state  errors  of  two  state  variables  are  known  for 
many  realizations  of  the  two  states,  then  a least  square  regression 
line  (see  Figure  (1))  can  be  described  by 

ZY  Y 

, ^1  ^^2 

Z— ' Y 7~ 

The  mean  values  of  the  errors  are  clearly  assumed  to  be  zero 
(i.e.,  no  biases  exist). 


FIGURE  1.  LEAST  SQUARES  REGRESSION  LINE. 
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For  a given  actual  value  of  Yui,  an  estimate  of  can 

be  made  from  the  equation  of  the  least  square  regression  line 
according  to  equation  (4) . The  "goodness"  of  the  estimate  depends 
on  how  well  Yu2  is  correlated  with  Yuj.  If  the  correlation  were 
perfect  (i.e.,  ri2  = then  there  would  be  no  scatter  in  the  data 
points  in  Figure  (1).  All  the  data  would  be  directly  on  the  linear 
regression  line.  Conversely,  a poor  correlation  would  be  seen  as 
much  scatter  of  the  data  about  the  regression  line  in  Figure  (1) , 
and  the  estimate  of  Yu2  from  the  regression  equation  would  not  be 
expected  to  be  very  good. 

Equation  (4)  can  be  rewritten  in  terms  of  the  covariance 
from  equations  (1)  and  (2)  as; 

est.  ^ 


A quantitative  measure  of  the  data  scatter  about  the 
regression  curve  is  defined  as  a "Standard  Error  of  Estimate" 


21 


(6) 


The  difference  of  the  two  terms  in  equation  (6)  is  the  difference 
between  each  data  point  and  the  regression  curve  approximating  the 
data.  The  standard  error  of  estimate  is  similar  to  a deviation  of 
the  data  about  the  regression  curve.  If  lines  are  constructed 
parallel  to  the  regression  line  of  Yu2  regressed  on  Yui  at  vertical 
distances  ^21,  2o2i,  etc.,  from  it,  then  for  large  N there  would  be 
included  between  these  lines  approximately  68%,  95%  etc-  of  the 
data  points.  Hence  ^^21,  is  analogous  to  the  standard  deviation  of 
normally  distributed  data.  The  Yu2  data  mean  for  any  given 
Yui  and  very  large  numbers  of  sample  data  lies  on  the 

regression  curve. 


A generalized  form  of  the  correlation  coefficient  given  in 
equation  (3)  is: 


Predicted  Variation 


yTotal  Variation 


and  for  linear  regression  it  can  be  shown  that 

? 


E[C“3)a  - ' -EC“2), 


(7) 


(8) 
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hence  with  (8) , (7)  and  (6) 

o = S (1  - r 

^22  ^21  ' 


B.  Multiple  Linear  Regression 

When  multiple  correlation  exist  among  a set  of  random 
variables  as  is  the  case  when  given  an  N x N covariance  matrix, 
multi-dimensional  regression  equations  can  be  constructed. 
Generalizing  equation  (5)  to  the  multi- dimension  equation: 




Where  the  an  coefficients  are  found  from  simultaneous  solution  of 
the  algebraic  equations: 


^1^11 


^1^21 


+ ^2^12 


^^2^22 


^3^13 


®3^23 


+ • • • + ®n-l^l(n-l) 

+ . . . + a . r , , . 

n-l^(n-l) 


♦ ^2’^  (n-1)  2 * (n-1)  3 • 


* (n-l)’^(n-l)(n-l)  '^(n-Dn 


and  the  standard  error  of  estimate  about  the  multiple  regression 
line  is  an  expanded  form  of  equation  (9) : 

% " ®nn  " ^n,l  ^1  " ^n,2  ^2  " * * ’ ^n,  (n-1)  ^n-l)^^%p) 

III.  INVERSE  REGRESSION 

A.  Method  Description 

If  it  is  desired  to  find  an  ensemble  of  state  vectors 
which  include  random  errors  and  if  it  is  also  necessary  that  the 
random  errors  included  in  the  state  vectors  have  the  proper 
statistical  cross  correlations,  then  it  may  be  possible  to 
construct  such  an  ensemble,  given  the  covariance  matrix  of  errors. 

The  process  would  be  an  inverse  regression  scheme  as  follows. 

Making  the  assumption  that  all  errors  are  normally 
distributed  about  linear  regression  lines  it  is  possible  to 
construct  a set  of  regression  lines  and  then  statistically 
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"scatter"  the  errors  about  the  line  using  a zero  mean  normally  ^ 
distributed  random  number  generator  with  a variance  equal  to  a 
(equation  12)  . The  errors  so  generated  would  be  added  to  the 
mean  value  of  the  state  vector  and  repeated  many  times  to  build 
an  ensemble  of  state  vectors  with  properly  correlated  errors. 

A consistent  scheme  of  selecting  the  order  of  regression 
should  be  used.  One  such  scheme  would  be  as  follows: 

1.  Select  the  two  elements  of  the  state  vector  that 
are  least  correlated  (i.e.,  smallest  ISnjtii  » N/M)  . 

2.  Continue  selecting  elements  in  a monotonically 
increasing  order  of  correlation. 


3.  Choose  the  selected  initial  two  elements  (from 
step  1)  and  construct  a regression  curve  from  equation  (5) 
similar  to  Figure  (1) . 

3a.  Find  the  so-called  actual  value  ^^'^1/ a from 
a zero  mean  normal  random  generator  (NRNG)  with  variance 


3b.  Solve  for  the  estimated  value  of  ^U2 

(y  \ 

4.  Introduce  scatter  to  V '^2ye  by  adding  an  error 
term  found  from  (NRNG)  using  variance  02  as  found  from  equation  (9) . 


Y 

Solve  for  the  so-called  actual  value  of  ^2 


(13) 


5.  Repeat  the  process  to  calculate  actual  values 
for  all  the  state  errors. 


Equation  (13)  is  a defined  value  of  a single  realization  of 
the  error  in  the  Y2  element  of  the  state  vector.  When  this  process 
is  repeated  for  each  element  of  the  state  vector,  it  is  necessary 
to  use  the  expanded  multiple  regression  equation  (10)  and  variance, 
equation  (12) . 


For  example  when  finding  the  third  error  term, 
from  (10) : 


(14) 


6 


I 


NSWC/WOL  TR  77-185 


where  (Yui)a  and  (Yu2)a  same  values  as  found  in  the 

precedinq  steps.  The  variance  used  to  introduce  "scatter"  into 
the  (Yu3jest.  of  equation  (14)  is: 


"3^  =[=33^''^  - 


®1  ■ ^32  ^2 


)‘i 


The  process  is  complete  when  all  N values  of  the  state  vector  have 
been  altered  by  an  error  term  as  in  equation  (la) . 

For  a given  covariance  matrix  of  errors,  the  Ajj  coefficients 
of  equation  (11)  need  only  be  calculated  once.  They  change 
only  when  the  covariance  matrix  changes.  However,  the  value  of 
the  "Ajj"  coefficient  are  different  for  each  value  of  N in 
equation  10.  A set  of  coefficients  correspond  to  each  value 
of  N which  is  the  order  of  the  given  covariance  matrix.  As 
many  realizations  as  desired  of  the  state  vector  can  be  generated 
for  a given  covariance  matrix.  Each  new  realization  is  constructed 
from  a new  (NRNG)  value  for  the  first  error  term  Yuj^. 

B.  Proof  of  Proper  Correlation 

It  remains  to  show  that  the  order  chosen  for  regressing  the 
variables  does  not  affect  the  statistical  cross  correlations  of 
the  errors  in  the  final  ensemble  of  state  vectors.  If  we  accept 
equations  (13)  as  a valid  method  of  introducing  "scatter"  or  random- 
ness to  regression  generated  errors,  then  it  is  sufficient  to  show 
that  cross  correlations  of  error  states  are  correct,  independent  of 
the  order  of  their  regression. 


('“n)a  - 


(NRNG) 


The  term  (NRNG)gj^  is  a single  realization  of  a normal  random  number 
generator  with  standard  deviation,  Oj.. 

Let  us  arbitrarily  choose  two  elements  of  an  error  state  vector 
and  Y^2  and  regress  Yu2  on 


N.  -fe-f-.) 


Where /y  \ is  chosen  fr 

(“ija 

02  being  found  from  (9) . 


is  chosen  from  (NRNG)  1/2 


Then 


(NRNG) 
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The  cross  correlation  between  ^ and(Yu2)a  ^ large  number 

of  realizations  of  equation  (18)  is: 

<“>12  “LC'-Oa  (S)a  /•> 

from  (17)  and  (18) 

<“>12  0''“l)a  * <™N<=>02  /'* 

- ^12  /" 

The  term  on  the  right  is  equal  to  zero  in  the  limit  since  each 
term  of  the  product  has  zero  mean  and  the  terms  are  uncorrelated. 
Hence  as  expected: 

<“>12  ■ ®12  =!](''“])  a (S)  a /«  <^ 


From  (10)  the  regression  equation  for  a third  variable  is: 
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This  quickly  reduces  to: 


. 1/2  „ 1/2  , . =33^^^  ®12  . (NBNG)„3 

- b33  A3  + - jyj  A2  + - 


’22 


(24) 


Since  the  last  term  is  zero  in  the  limit  for  a product  of 
uncorrelated  zero  mean  variables,  (24)  is: 


N 


S 

a _ „ 1/2  „ 1/2  ^33  ^12  , 

- ■ ®11  ®33  ^1  ^ _ T75  ^2 

^22 


(25) 


The  A3  arid  A2  coefficients  are  found  from  (10)  to  be; 


^13  " ^12  ^23  , ^23  " ^12  ^13 

^1  = — ; 2 — ' ^2  = ^5 — 


1 - r 


12 


1 - r 


12 


(26) 


Substituting  (3)  and  (26)  into  (25) ; 


Hence,  even  though  YU3  was  the  third  variable  chosen  it  is 
properly  correlated  with  the  first  variable.  In  a similar  fashion, 
the  third  variable  can  be  shown  to  be  properly  correlated  with  the 
second  variable.  This  process  can  be  continued  for  all  subsequent 
elements  of  the  state  vector  regardless  of  their  order  of 
regression.  Therefore,  the  order  chosen  to  regress  the  state 
vector  errors  is  arbitrary  and  has  no  adverse  effect  on  the 
cross  correlations  of  the  ensemble  of  state  vectors  so  generated. 


IV.  CODE  DESCRIPTION 

The  inverse  regression  method  described  has  been  coded  for 
use  in  a trajectory  program.  If  a state  error  covariance  matrix  is 
available  at  some  point  in  a vehicle's  trajectory  and  it  is 
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1 


desired  to  complete  the  trajectory  with  a Monte-Carlo  error 
analysis,  then  it  is  necessary  to  determine  an  ensemble  of  initial 
state  vectors.  One  state  vector  from  the  ensemble  is  used  as  an 
initial  condition  for  each  trajectory  of  the  Monte-Carlo 
series. 

A listing  is  enclosed  as  Appendix  B.  for  the  subroutines 
that  read  in  the  covariance  matrix  and  then  determine  the 
coefficients  and  standard  deviation  required  to  solve  for  state 
error  vectors.  Subroutine  COVSET(M)  stores  in  common  the  "A" 
coefficients  and  "S"  variances  of  equation  (10) . Also  stored  are 
the  standard  deviations  of  equation  (12)  for  the  M-dimensional 
square  covariance  matrix.  These  stored  parameters  along  with 
a normal  random  number  generator  are  sufficient  to  form  a state 
error  vector  by  repeatedly  solving  equation  (10)  for  M values  of 
the  state  est.  This  error  state  determines  the  initial  condition 
for  one  trajectory  realization.  Additional  realizations  of  the 
Monte-Carlo  series  can  be  found  with  the  same  stored  coefficients 
(i.e.,  these  need  not  be  calculated  again).  The  normal  random 
number  generator  with  previously  calculated  standard  deviation 
supplies  the  required  variation  of  the  state  vector  for  subsequent 
trajectories . 

Subroutine  FILL  uses  a system  supplied  function  (MAM)  to 
solve  the  set  of  simultaneous  equations  (11) . Any  simultaneous 
algebraic  equation  solver  may  be  substituted  here.  Subroutine 
IMOD  actually  generates  the  desired  error  vectors. 


V.  CONCLUSIONS 

A method  has  been  described  which  allows  mixed  modes  of 
error  modeling  for  a single  event.  The  method  defines  the 
interface  between  a covariance  propagation  error  model  and  a 
Monte-Carlo  deterministic  error  model.  Given  an  arrival 
uncertainty  covariance  matrix  for  a physical  system  at  some  point 
in  time,  it  is  possible  to  construct  an  ensemble  of  state 
error  vectors  which  may  be  used  as  initial  conditions  for  the 
error  state  in  order  to  exercise  a Monte -Carlo  model  of  the  system. 
The  error  states  so  generated  have  been  shown  to  be  properly 
cross  correlated  with  each  other  and  have  the  correct  distribution 
of  magnitudes.  A FORTRAN  computer  program  has  been  written 
to  facilitate  implementation  of  the  method  within  trajectory  codes. 
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APPENDIX  A 


SYMBOLS 


a 

cc 

N 

NRNG 

r 


®NN 


Regression  coefficients 
Error  covariance 

Number  of  samples  or  summation  index 
Zero  mean  normal  random  number  generator 
correlation  coefficients 
Variance  of  error 

Covariance  of  errors 

Element  of  state  vector 

Known  part  of  state  vector  element 

Unknown  part  of  state  vector  element 


o 

E 

Subscripts 
est.  or  e 
actual  or  a 


Standard  deviation  for  gaussian  distribution 
Summation  symbol 

Regression  estimated  value  of  an  error 
Determined  actual  value  of  an  error 


A-l/A-2 


d d ood  n n n ono 
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APPENDIX  B 

COMPUTER  CODE  LISTING 
SI'HROUTINF  ^n^/‘;PT(^) 

HPiiNsvcLn  '•'v  y? 

FOP  * GIVPM  ^•^NT^ pt  '■»  S*^T»  THf  fOLLO'^ING  I MF  opt^  A T 1 0^‘ 

NPEH  «E  CAl  rci  'JTt-i-  OM.Y  O^-Cf  . 

CALCULATE  CFPTMM  COF  c r IC  I rh.JS  FOOM  ftt^KlVAl  UNCFPTATNTY 
COVAPTANCF  mATpIX  (T^'TTI^L  rf>Nl- 1 T I ON«,  K Ow  EACH  P^l^) 

storage  IM  THF  < tPFAY  SECpc  AT  X(l)  ANO  TS  ENOEH 

storage  LOCATION'S  = *A*  STOwAGF  ♦ SlGt-o  STORAGE  ♦ S(N)  STOOAGE 
= ( (R-1  ) ♦ (R-?) 

X ARRAY  NO'm  S"RFICIFNJT  T(j  HAMOLE  ' 1 >>  X 1? 


COMMO^•  X(RA) 


OIMFnJFION  .CMA*  (■»,?)  .A*<HAY(P.?)*IC(?) 

?00  format (5F1S, 7. 5X) 

?01  FORMATdH  ^SFl^i.T.RX) 

?0?  FORRAT(lHl) 

X (?)=FLOAT (R) 

READ  IN  covariance  RATPIX 
READ (S, POO) SMM 

WHITE (ft, 20?) 

WRITE (ft, 201 ) SMM 

TO  SIMPLIFY  C0MP(JTATI0N5»  niA(iONAL  ELF^tNIS  Oh  SRM  MATRIX  ARE 
SORT  oF  original  COVAPTAMrF  OIAGONaL  ELF^EnTS. 

CALCULATE  CORPPLATION  COfFP  ICIFNT^  ^ATmIX  FRO-"  ARRIVAI  COVARIANCE 
on  10  K=1,M 
no  10  J=1.M 
IF  (K.GE.J)('0  TO  10 

CALClILAIF  the  *o«  VAI  tiPS  (CORHtLATIOf  COEFFICIENTS) 

correlation;  C0FFFTC1F^:T  .XiT'^IX  RFCOVES  AS  V^E  yRTTL  ON  TOP  OF 

original  COVARIANCE'  CATRIY  TO  SAVE  STfPAGK. 

SRM  ( J,^)  =S'^M  ) / (SRR  ( J»  J)  (K,K  ) ) 

SmM(K,J)=SMN'(J,K) 

10  CONTINUE 

DO  11  J=1,M 

PUT  THF  S(N)  values  IN'  FTr.pAOE,  THESE  aE£  ThE  OIAfONAL  ELF^tHTS  OF 
THE  COVAHIAMCF  MATRJX  FOItM  TO  THE  STANDARO  DEVIATION,  FIRST 
LOCATION  IS  Y (?) . 

X (?*J)=SMM(J, J) 

SET  THF  OIAGOMAL  flFmF'MTS  OF  THF  CORRELATION  vaTPIY  TO  l.C 

11  SNM(J,J)=1, 

X(2)=M  ANO  X ( 1 ) =S  ( 1 ,?) /s  ( 1 , ] ) , ARE  roE  PFH'ieNTS  FOR  Yt«?  FFTIMATF 
X(l)=SMM(l.?)*x (4)/X(3) 

Nil  i;  R3  = Mr  = M + 2 1. 


FIRST 


CALCULATE  THF  *A*S  AN'O  PUT  TMTil  STORAr-E,  A(l)  ANU  A(2)  FIRST. 
these  are  THE  UNKN'Orn  VARIARLFS  of  The  SIRUL  TANFOUS  ALGFHRAIC 
FOUATIOMS,  SUHRf)UTlN'F  FILL  SOI  VES  THFSF  FOI|ATIun!«.. 

BMM=SOLUTTONf  VFCTCiRS  = «^a*S,  , FVM=CuRRfcLAT  1 On  COFFUCIEMT  maTPIX, 
MsOROFR  of  Smm,  m=nUnhfp  of  siMUl  TanEOUS  FOIIaTioniS  To  HF  SOLVE!' 

IN  THIS  PASS  through  FILL.  ARRAY  ANO  rC=l.UMMY  •MATRIX  THAT  MUST  BE 
OTMENSIONFO  TM  roVSET.  OI'^EnSIOn  TO  m-1 


BtSIIAVAIlABlE  COPY 


BEST.  AYAIlABiE  COPY 
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I 


C 


c 

c 

c 

r. 


c 


c 

c 

c 


c 

c 

c 

c 

r 


c 


c 

c 


c 


r 


n 


14 


21 


CALL  FILL  ir  ) 

on  13  ILOs],M 

MMrlLn^Mfl 

STORE  SOLUTIOM  VFfTnpc;  FILL*  FI^'ST  LOCATIOM  ic  , 

X (MM) =HMM ( ILO) 

CONTlAiOE 

IF (N.EO.Ml ) GO  TO  14 
MH  = MM 


60  TO  1? 

calculate  m-1  PFRMANFkT  VAlUIS  OF  mof-f  , SIGMa(?)  !«;  FIRST  aI 

LOCATION  X (MM+l ) =x (3  + M* (M-) ) ♦ ( i-p) 

calculate  <:V,MA(;^) 

X {MM*l ) =X (4) * ( 1 .-Smm ( 1 ,?)*«?) *«.s 
calculate  sigma  (3)  TO  SIG**A  ("•) 

DO  21  K=?*M1 
RFST=0. 

DO  22  L=1*K 

REST  = SMM  (K  + i ,1.)  «X  (MC*L  ) ♦RFST 

CONTINUE 

MC=MC+K 


PUT  CALCULATFO  SIGMA  s INTO  X 
X(MM  + k)=X(K*3)*(1  ,-PtST 


CONTINUE 

WITH  sigma  S*  *A«S.  ami.)  SCO  S»  fr«0>'  veCluRS  VAY  af  f;KMFP«T*-f' 
WITH  normal  PANOOM  M(G‘«MFH  oFM^RATOP.  A SL'Hwoi  ;T  I N'  ^ 0’^  THic 


called 
RE 'I'rn 


"IMOD"  IS  INCLUOFO, 


END 


SURROUTINF  F I LL ( « . PMM .M , N , A Foa Y . IC ) 

BRUNSVOLD  may  77 

THIS  SURROUTINF  FINDS  ThF.  fql  I iT  IO-"  VECTOR  FOw  ThF  COt  F F T 0 1 F 'Tc 


DIMENSION  PMM(M,v) . AHway (M«M) .H (N) , IC (a ) 

DO  IS  K=1.N 

GENERATE  RIGHT  SIDE  OF  ALGFRPaIC  EnUATIONS  F PO^'  CC'^f  1 1 A T T 0^.  mATfIy 
B (K)  =PMM  (K,N-»l  ) 

DO  IS  J=1,N 

GENERATE  COEEEICIENTS  OF  AlGFPPATC  FOUaticivS  f pom  rOPt-tlATION 
MATRIX  AND  PUT  INTO  AHRAY. 

APRAY(J»K)=PMM(J,K) 

APRAY (K tJ) =APPAY (J.K) 

15  CONTINUE 

mam  IS  system  SIJHPOUTTNF  That  SOLVFC  I InEak  SF  T of  ALGFHPAIC  funs, 
call  mam ( ADPAY*N.N*P* 1 » IC* 10) 

IF(ID.EQ.2)G0  TO  v>0 
GO  TO  16 
20  wPiTE(6*inn) 

ion  FORMAT  (IN  »2SHFAIIUPE:  TU  «UF<PO"TInE  FUL) 

STOP 

RFTUPM  TO  rOVSFT  WITH  foLUTIOF'  VFf  tor  t ■ . 

16  RFTURN  $ FND 


I 
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SDBPnUTINF  IMHO 

r this  tnitial  vont/t  e solvi^s  fop  thf  initial  fppops  tw  tme  ptatf  vff 

COMMON  X(PO) 

DIMENSION  Y{Q) 

M=IFIX (X (P) ) T nA=n. 

C FIND  THF  PBQPPP  INDEX  IN  Thf  > fiopny 

J=0  S Ml=M-l 
DO  S I=?.M 
J=I^J 

S COMTIMIJE 

GENERATE  A NORMALLY  OISTRIBUTFO  PANDON  NJ^hfR*  (R^;)  * OF 
SIGMA=1.  AND  MFAN=0. 

?S  CALL  RANNIIM(OA.1,*0.»P»PN) 

GENERATE  THE  FIRST  STATE  PCROR  Yd). 

Y ( 1 ) =PN*X (3) 

CALL  PAMNI)M(OA,1.,0,.P,RN)  [ 

GENERATE  THF  SFgOh’O  STATE  frR(p,  | 

Y(?)=X(l)*Yd)+  n»X(n'M+1)  I 

mc=m+? 

SOLVE  FOR  The  PFmaINImG  Y ERRORS 
DO  in  K=?*Mi 
REST  =0. 
nr>  ?0  L = 1»K' 

RFST= (X (HC+L )*Y(L)/X(L*?)) ♦REST 
?n  CONTIA'UF 

call  panNIIm(0a.1  ,,o,.r,PN) 

Y(K^1)=X(K*3)  »PKST  + R|vI*X  (M^  ♦k  ) 

Mr=MC*K 
10  CONTINUE 

THE  Y VFCTOh  PF^lFRATFP•  IS  O^iF  frrnr  V^OTOP  pFALI7ATT0n  TO  BE  AOnFP 
TO  means.  AOOITinMAL  REAL  I RATIONS  CAf  HE  GENERATED  WITH  OTHFP  SET 
OF  PA‘'I(10M  MIIM°FPS*  RM'>S. 
k<PITE  (p*inn)  Y 
inn  format ( iHn.PFi o.s) 

PFTUR»'  S END 
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Attn:  Engineering  Libary 

Dept,  of  Aeronautics  and  Astronautics 
Seattle,  Washington  98105 

West  Virginia  University 
Attn:  Library 

Morgantown,  West  Virginia  26506 

Federal  Reports  Center 
University  of  Wisconsin 
Attn:  S.  Reilly 

Mechanical  Engineering  Building 
Madison,  Wisconsin  53706 

Los  Alamos  Scientific  Laboratory 
Attn:  Reports  Library 

P.  O.  Box  1663 

Los  Alamos,  New  Mexico  87544 

Institute  for  Defense  Analyses 
Attn:  Classified  Library 

400  Army-Navy  Drive 
Arlington,  Virginia  22202 
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University  of  Florida 

j Attn;  Dept.  Engineering  Science 

I Dr.  B.  M.  Leadon 

] • Gainesville,  Florida  32601 

I Keunan  Sciences  Corporation 

. Attn:  Library 

Mr.  D.  Foxwell 
P.  0.  Box  7463 

j Colorado  Springs , Colorado  80933 

t 

, Rockwell  International  Corporation 

Technical  Information  Center 

“ 4300  E.  Fifth  Avenue 

f Columbus,  Ohio  43216 

i M.I.T.  Lincoln  Laboratory 

I Attn;  Library  A-082 

Dr.  A.  B.  Wardlaw 
P.  O.  Box  73 

Lexington,  Massachusetts 

i The  RAND  Corporation 

! 1700  Main  Street 

Library  - D. 

Santa  Monica,  California  90406 

Aerojet  Electrosystems  Co. 

Engineering  Library 
1100  W.  Hollyvale  Avenue 
Azusa,  California  91702 

The  Boeing  Company 
Attn:  87-67 

P.  O.  Box  3999 
Seattle,  Washington  98124 

United  Aircraft  Corporation 
Attn;  Library 
400  Main  Street 

East  Hartford,  Connecticut  06108 

Hughes  Aircraft  Company 
Attn;  Company  Tech.  Doc.  Center 
6/Ell,  B.  W.  Campbell 
Centinela  at  Teale 
Culver  City,  California  90230 
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t 
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Copies 

Lockheed  Missiles  and  Space  Co.,  Inc. 

Attn:  C.  Grolemund 

P.  O.  Box  504 

Sunnyvale,  California  94086 

Lockheed  Missiles  and  Space  Co.,  Inc. 

Attn:  Technical  Information  Center 

3251  Hanover  Street 
Palo  Alto,  Califonria  94304 

Lockheed-California  Co. 

Attn:  Central  Library  Dept. 

Burbank,  California  91503 

Vice  President  and  Chief  Scientist 
Dept.  03-10 

Lockheed  Aircraft  Corp. 

P.  O.  Box  551 
Burbank,  California  91503 

Martin-Marietta  Corporation 
Attn:  Science-Technology  Library 

(Mail  No.  398) 

P.  0.  Box  988 
Baltimore,  Maryland  21203 

Martin  Marietta  Corp. 

Orlando  Division 
P.  0.  Box  5837 
Attn:  Mr.  H.  J.  Diabolt 

Orlando,  Florida  32805 

General  Dynamics 
Attn:  Research  Library  2246 

P . 0 . Box  748 
Fort  Worth,  Texas  76101 

Calspan  Corporation 
Attn:  Library 

4455  Genesee  Street 
Buffalo,  New  York  14221 

Air  Force  University  Library 
(SE)  63-578 

Maxwell  Air  Force  Base  • f 

Alabama  36112  j 
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McDonnell  Douglas  Corp. 

Attn:  R.  D.  Detrich,  Dept.  209  Bid.  33 

P.  O.  Box  516 

St.  Louis,  Missouri  63166 

McDonnell  Douglas  Astronautics 
Company  - West 

* Attn;  J.  S.  Murphy,  A3 

368,  B4A0 
5301  Bolsa  Avenue 

Huntington  Beach,  California  92647 

Fairchild  Hiller 
Engineering  Library 
Republic  Aviation  Division 
Farmingdale,  Newyorh  11735 

General  Applied  Sciences  Labs,  Inc. 

Attn:  Dr.  F.  Lane 

L.  M.  Nucci 

Merrick  and  Stewart  Avenues 
Westbury,  Long  Island 
New  York  11590 

General  Electric  Company 
R&D  Labs  (Comb.  Bid.) 

Attn:  Dr.  H.  T.  Nagamtsu 

Schenectady,  New  York  12301 

The  Whitney  Library 
General  Electric  Research  and 
Development  Center 
Attn;  M.  F.  Orr , Manager 
The  Knolls,  K-1 
P.  O.  Box  8 

Schenectady,  New  York  12301 

General  Electric  Company 
Missile  and  Space  Division 
Attn;  MSD  Library 
P.  O.  Box  8555 

Philadelphia,  Pennsylvania  19101 

» General  Electric  Company 

AEG  Technical  Information  Center,  N-32 
Cincinnati,  Ohio  45215 
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General  Electric  Company 

Reentry  & Environmental  Systems  Division 
Attn:  J.  Immel 
3198  Chestnut  Street 
Philadelphia,  Pennsylvania  19101 

AVCO- Everett  Research  Lab. 

Attn:  Library 

2385  Revere  Beach  Pkwy. 

Everett,  Massachusetts  02149 

LTV  Aerospace  Corporation 
Vought  Systems  Division 
Attn:  Unit  2-51131  (Library) 

P.  O.  Box  5907 
Dallas,  Texas  75222 

LTV  Corporation 
Attnj  MDS-T-Library 
P.  0.  Box  5907 
Dallas,  Texas  75222 


Northrop  Norair 
Attn:  Tech.  Info.  3343-32 

3901  West  Broadway 
Hawthorne,  Califonria  90250 

Government  Documents 
The  Foundren  Library 
Rice  Institute 
P.  O.  Box  1892 
Houston,  Texas  77001 

Grumman  Aerospace  Corporation 
Attn:  Dr.  R.  E.  Melnik 

Bethpage,  Long  Island 
New  York  11714 

Marquardt  Compemy 
Attn:  Library 

P.  0.  Box  2Ul3 

Van  Nuys,  California  91409 

ARDE  Associates 

Attn:  Librarian 

P.  O.  Box  286 

580  Winters  Avenue 

Paramus,  New  Jersey  07652 
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Aerophysics  Company 
Attn:  Mr.  G.  D.  Boehler 

3500  Connecticut  Ave. , N.W. 

Washington,  D.C.  20003 

Aeronautical  Research  Associates  of  Princeton 
Attn;  Dr.  C.  duP.  Donadson 
50  Washington  Road 
Princeton,  New  Jersey  08540 

General  Research  Corporation 
Attn;  Tech.  Info.  Office 
5383  Hollister  Avenue 
P.  O.  Box  3587 

Santa  Barbara,  California  93105 

Sandia  Laboratories 
Attn:  Dr.  C.  Peterson 
Dr.  G.  W.  Stone 
Dr.  K.  Touryan 
Dr.  R.  Eaton 
Kurt  Putz,  Div.  1333 
Box  5800 

Albuquerque,  New  Mexico  87115 

Hercules  Incorporated 
Attn:  Library 

Allegany  Ballistics  Lab. 

P.  O.  Box  210 

Cumberland,  Maryland  21502 

General  Electric  Company 
Attn;  Dvae  Hovis,  Rm.  4109 
P.  O.  Box  2500 

Daytona  Beach,  Florida  32015 
TRW  Incoroorated 

Attn;  Tech.  Library/Doc.  Acquisitions 
1 Space  Park 

Redondo  Beach,  California  913278 

Stanford  Research  Institute 
Attn;  Dr.  G.  Abrahamson 
333  Ravenswood  Avenue 
Menlo  Park,  California  94025 

Hughes  Aircraft  Company 
Attn:  Tech.  Library,  600-C222 
P.  O.  Box  3310 

Fullerton,  California  92634 


15 


NSWC/WOL  TR  77-183 


Copies 


Westinghouse  Electric  Corp. 
Astronuclear  Laboratory 
Attn:  Library 
P.  0.  Box  10864 

Pittsburgh,  Pennsylvania  15236 

University  of  Tennessee 
Space  Institute 
Attn:  Prof.  J.  M.  Wu 

P.  0.  Box  10864 

Pittsburgh,  Pennsylvania  15236 

CONVAIR  Division  of  General  Dynamics 

Library  and  Information 

P.  O.  Box  12009 

San  Diego,  California  92112 

CONVAIR  Division  of  General  Dynamics 
Attn:  Research  Library 

P.  O.  Box  80986 

San  Diego,  California  92138 

AVCO  Missiles  Systems  Division 
Attn:  E.  E.  H.  Schurmann 

N.  Tyson 
H . Rosenbaum 
201  Lowell  Street 
Wilmington,  Massachusetts  01887 

Chrysler  Corporation 
Space  Division 
Attn:  N.  D.  Kemp  2910 

E.  A.  Rawls,  Dept.  2920 
P.  0.  Box  29200 

New  Orleans,  Louisiana  70129 

General  Dynamics 
Pomona  Division 

Attn:  Tech.  Doc.  Center  (6-20) 

P.  O.  Box  2507 

Pomona,  California  91766 

Ford  Aerospace  & Communications  Corp. 
Ford  & Jcunboree  Roads 
Attn:  Dr.  A.  Demetriades 

Newport  Beach,  California  92660 
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Raytheon  Company 
Attn:  D.  P.  Forsmo 

Missile  Systems  Division 
Hartwell  Road 
Bedford,  Maine  01730 

JCBW  Defense  & Space  Systems  Group 
Attn:  M.  W.  Sweeney,  Jr. 

Space  Park  Drive 
Houston,  Texas  77058 

Marine  Bioscience  Laboratory 
Attn:  Dr.  A.  C.  Charters 

513  Sydnor  Street 
Ridgecrest,  California  93555 

Applied  Mechanics  Review 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78228 

American  Institute  of  Aeronautics  and 
Astronautics 
Attn:  J.  Newbauer 

New  York,  New  York  10019 

Technical  Information  Services 
AIAA 

Attn:  Miss  P.  Marshall 

750  Third  Avenue 

New  York,  New  York  10017 

Faculty  of  Aeronautical  Systems 
University  of  West  Florida 
Attn:  Dr.  R.  Fledderman 
Pensacola,  Florida  32504 

Space  Research  Corporation 
Chittenden  Bank  Building 
Attn:  Library,  J.  A.  Finkel 

North  Troy,  Vermont  05859 

The  Aerospace  Corporation 
Attn:  E.  Ndefo 

J.  M.  Lyons 
P.  O.  Box  92957 

Los  Angeles,  California  90009 
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Notre  Dame  University 
Department  of  Aerospace  Engineering 
College  of  Engineering 
Library 

Notre  Dcime,  Indiana  46556 

Acurex  Corp.  Aerotherm 

485  Clyde  Avenue 

Attn:  M.  Abbett 

Mt.  View,  California  94042 

Mathematics  Research  Center 
Attn:  Prof.  Seymour  V.  Par ter 

U . S . Army 

University  of  Wisconsin 
Madison,  Wisconsin  53706 

TRW  Defense  & Space  Systems  Group 
One  Space  Park 
Attn:  B.  Pearce 

Aerodynamics  Dept. 

Garth  W.  Lippmann,  Bldg.  R-5,  Rm.  2230 
Redondo  Beach,  California  92078 

Lockheed  Missiles  and  Space  Co.» 

Continental  Bldg.,  Suite  445 

Attn:  R.  Fortune 

El  Segundo,  California  90245 


